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I. INTRODUCTION 



Among other uses, the heat kernel [l[ is a tool to deal with one- loop effective actions in quantum field theory. 
The effective action, the trace of the logarithm of the fluctuation operator Q, suffers from ultraviolet divergences, as 
well as many-valuation and anomalies. As noted in |3j the heat kernel has the virtue of being one- valued, free from 
ultraviolet divergences and gauge invariant. 

The heat kernel finds a number of applications: study of spectral densities of Klein-Gordon operators, proof of 
index theorems [1, [E| , to compute the (- function Q and the anomalies of Dirac operators Q , to deal with chiral gauge 
theories Q and models of QCD Q, to the Casimir effect [l(|, to compute black hole entropies [ll|, etc. 

Except in very particular manifolds, the heat kernel is expressed by means of asymptotic expansions. The Seeley- 
DeWitt expansion [12, HH, is in powers of the proper time, and is available to rather high orders in several setups, 
including curved spaces with and without boundary, and in presence of non Abelian gauge fields and non Abelian 
scalar fields, using different methods fll. [I3 - l20l | . 

To study quantum field theory at finite temperature one can use the imaginary time formalism, with compactified 
Euclidean time [U [22[ . This introduces a modification in the heat kernel coefficients. Early attempts to compute 
those coefficients were made in [23l |24|. However, ad hoc assumptions made in those calculations (essentially what 
we call the quenched approximation below) lead to expressions in conflict with explicit results derived for particular 
settings [25ll26j. The first systematic, and fully gauge covariant, calculation of the heat kernel at finite temperature 
was presented in [2?], H^|. There it was found that, besides the usual covariant derivatives, the Polyakov loop, Q(x), 
was also present in the expressions (consistently with (2f| [26j]). This is to be expected, since the Polyakov loop is the 
other natural gauge covariant construction allowed at finite temperature. This is not just a technical nicety, in fact, 
nowadays the gluonic Polyakov loop in QCD at finite temperature plays a prominent role as a relevant order parameter 
of confinement in the very successful Polyakov-Nambu-Jona-Lasinio models [29l - [3l| . The Polyakov loop appears 
automatically in any gauge covariant computation at finite temperature, and solves long standing paradoxes related 
to gauge invariance due to naive perturbative expansions [32h341 |. Moreover it is the only way a chemical potential 
could appear in the effective action. Indeed, the chemical potential is obtained by the shift Aq(x) — > Aq(x) — fi, where 
H is a constant real c-number. This has no effect in [Do, ], but it shows through the Polyakov loop dependence due 
to fi(x) ->■ e^n{x) pi [H. 

The results of (27|,[28| refer to the usual heat kernel expansion. That is, the coefficients are classified according to the 
dimension of the operators they carry (this classification holds at zero or finite temperature, and at zero temperature 
is equivalent to an expansion in the powers of the proper time). In [35] an expansion of the (zero temperature) heat 
kernel based on the number of covariant derivatives was carried out. This is a resummation of the usual expansion 
in which each coefficient has a fixed number of covariant derivatives but any number of scalar fields. The extension 
to curved space-time was made in [36j . In the present work we compute, for the first time, the heat kernel at finite 
temperature within the covariant derivative expansion. 

The results for the heat kernel at finite temperature of [27l |28| were obtained using a rather cumbersome method. 
Essentially it was a mixture of (already known) zero temperature coefficients for the spatial covariant derivatives plus 
the method of symbols [37l |38| for the covariant time derivative. In this approach some work is required to bring 
the expression to a manifestly gauge covariant form, involving the Polyakov loop. This is largely improved in the 
present paper. The new idea presented here is based on extending the method of covariant symbols, introduced by 
Pletnev and Banin 39], to the finite temperature case. The original method was devised for zero temperature, and 
so it assumed a continuous frequency variable. We adapt here the method so that it applies also for the discrete 
Matsubara frequencies. The Polyakov loop is accommodated in a natural way in the new approach. By means of this 
new technique, the calculation of the heat kernel at finite temperature, or other quantities like the effective action, can 
be done with manifest gauge covariance at each step. The method applies to general pseudo-differential operators. 
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In loop momentum integrals, the spatial components are continuous, but the frequency becomes discrete as a 
consequence of periodicity. This is equivalent to introducing a weight function in momentum space which consists of 
a family of Dirac deltas with support at the Matsubara frequencies. Here we find the remarkable result that much 
of the formalism goes through also for completely general weight functions h(p) in momentum space. This allows to 
obtain expressions which look Lorentz covariant (prior to momentum integration). The finite temperature case can 
be obtained from the generic one by replacing h(p ) by its Matsubara version. As a third contribution of this work, we 
adapt Chan's formula for the effective action [40| to such /i-spaces, and so in particular to finite temperature. (This 
automatically implies the corresponding result for the heat kernel.) The existence of Chan's form in such general 
setting is far from obvious a priori since the original construction by Chan relied heavily on integration by parts and 
averages in momentum space. These tools are not available in the presence of a generic weight h(p). 

The paper is organized as follows. In Scction|TT]we make a summary of previous results and techniques and develop 
the new method of covariant symbols valid at finite temperature, eq. (|2.40[) . In Section [IIII we present explicit results 
for the strict covariant derivative expansion of the heat kernel at finite temperature, to third order for the diagonal 
matrix elements, eq. p.26|) . and to fourth order for the trace, eqs. (|3.31j) and ()3.32p . Non stationary and non Abelian 
configurations are assumed throughout. In Section llVI we extend the gauge covariant technique to /i-spaces and use 
it to obtain the very compact Chan's form of the effective action, eqs. (|4.22|) and (|4.23|h In Section [V] we summarize 
our conclusions. Some auxiliary material and results are given in the appendices. 



II. METHOD OF SYMBOLS 



A. General considerations 



Let us consider a theory of scalar fields in rf-dimcnsional Euclidean fiat space-time, coupled to external fields, 
including gauge fields. Typically 

C(x) = -(f>(x^K4>(x), K = D 2 + X(x), D li = d lt + A ll (x). (2.1) 

The external fields X(x) and A^ix) are matrices in internal space in general. For concreteness we assume that <j>(x) 
transforms in the fundamental representation of the gauge group, 4>(x) — > U~ l (x)4>(x). 
The corresponding partition function and effective action are 

Z = J V<jJVcf)e-S ddxC{x) = e~ r , r = Trlogif. (2.2) 

The effective action T is a functional of the external fields and diagrammatically Tr log corresponds to adding one-loop 
graphs with the field 4> running in the loop and any number of external legs attached to it. 

The operation Tr can be expressed as a trace on a single-particle Hilbert space where K acts. This Hilbert space 
includes space-time and also internal degrees of freedom: 



r = J d a xtv{x\\ogK\x), (2.3) 

|a;) is a basis of the space-time sector, 

(x\x'} = S(x — x'), x^,\x) = x^x), (2.4) 

and tr refers to the internal degrees of freedom. Likewise, under a variation of the gauge fields and the scalar field, 
one obtains the current and density, 

5T = J d d xtr{J li {x)5A li {x)+V{x)5X{x)), 
J»{x) = {x\{K-\D^\x), V{x) = {x\K- 1 \x). (2.5) 

These examples, as well as the heat kernel, exp(r_R"), to be considered later, illustrate the need for computing 
diagonal matrix elements of pseudo-differential operators. Taking coincident points amounts to integrate over the 
momentum of the loop. 

In view of the above, we consider a generic pseudo-differential operator 



f = f(D,X), 



(2.6) 
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constructed with the covariant derivative and other fields, X{x). These external field are bosonic. The quantum 
field running in the loop may be bosonic or fermionic. Under a gauge transformation, — > U~ 1 D^U , X — > U~ l XU , 
the diagonal matrix elements transform covariantly, (x\f(D,X)\x) —> U~ 1 (x)(x\f(D,X)\x)U(x). 

Our goal is to address the computation of the diagonal matrix elements of the pseudo-differential operator, 
(x\f(D, X)\x), and its trace, in a gauge covariant setting valid at zero or finite temperature. 



1. Covariant expansions at zero temperature 



In general the diagonal matrix element cannot be expressed in closed form. At zero temperature, a typical expansion 
to be applied is one based in powers of and of X(x). This produces an expansion in terms of local gauge covariant 
operators 

(x\f(D,X)\x)=Y,9xOx(x). (2.7) 

A 

Here 0\(x) includes all possible local gauge covariant operators constructed with and X. That is, with X, with 
the field strength tensor, 

F la , = [D li ,D v ], (2.8) 

and with their covariant derivatives. The coupling constants g\ depend on the concrete operator /. Often the terms 
are organized by dimensional counting in subsets of operators with a common dimension. An example is the standard 
heat kernel expansion 

- j^j, (l + rX + fa* + \X„ + i**) +•••). (2.9) 
We indicate covariant derivatives using the convention 1 

YlJ,llJ,2— Mra = [D/H ) ^Al2—Mn] ' (2-10) 

for any operator Yj with a (possibly empty) ordered set of Lorentz indices /. So, for instance, F a ^ u — [Z3 Q ,F M ^] and 
X a p — [D a , [Dp,X]]. 

Another expansion, which is the subject of this work, is the covariant derivative expansion, which is a resummation 
of the previous one: at a given order the number of is fixed while there can be any number of X. For Abelian X 
this is just of the form 



(x\f(D,X)\x) = ]T/ A (A(.*))0 A (x), (2.11) 

A 

where now Oa(^) contains only X with derivatives, and f\(X(x)) is a generic function of X. In the more general 
case of non Abelian fields, one can still express the expansion by means of labeled operators [35ll4l|: 

{x\f{D,X)\x)=Y,h{Xi{x),...,X n {x))O x {x). (2.12) 

A 

The idea is that 0\(x) is the product of n — 1 local covariant blocks and the i-th copy of A, Xi, is meant to act 
between the (i — l)-th and the i-th block. For instance 

J* e «i£,e<«-*>* dt = J* e^-V*' dtF%, = eSX x ~_ e x 2 2 F l»- ( 2 - 13 ) 

Here X\ is X acting at the left of F® and X2 is X acting at the right. Note that the labeled operators Xi can be 
treated as c-numbers since A1A2 = X2X1. 



Here and elsewhere Y denotes a generic operator. 



As an example, all the terms of the heat kernel with precisely one X MM can be collected in the form [35[ 
Expanding in powers of Xi gives back the standard heat-kernel expansion: 

e rX 1+e rX 2 2e rX 1 _ e rX^ ,2 3 X 

X^ = + To V 1 + A 2j H -A 



(X!-X 2 ) 2 T(Jd-X 2 ) 3 y W_ V6 12 v ; / MA1 

r 2 t 3 



-X^ + -{X,X^} + -.. . (2.15) 



3. Covariant expansions at finite temperature 

At finite temperature the space-time is R d_1 x S 1 , within the imaginary time formalism [2ll |22| |. The quantum 
field may be bosonic or fermionic, being respectively periodic or antiperiodic in time with period j3 = where T 
is the temperature. The external fields A^{x) and X(x) are bosonic and hence periodic. The gauge transformations, 
U(x), are also periodic. 

The expansions in eq. (|2.7I) and (I2.12[) refer to zero temperature and they have to be modified at finite temperature. 
In fact, at finite temperature there are two gauge covariant constructions with the operator D$, namely, the covariant 
derivative [D , } and the Polyakov loop, 

n{ x ) = Pe -LT" A ^ dt . (2.i6) 

The Polyakov loop here is not traced, it is a matrix in internal space and P refers to path ordered product. Also, the 
integral starts at Xq rather than zero. The Polyakov loop so defined is gauge covariant at x: 

Q.(x) -> U-\x)n(x)U(x). (2.17) 

Q(x) is also periodic in xq. In practical terms, VL(x) behaves as a local field. This operator appears through Dq due 
to the relation [27| 

e -P Do = Q(x). (2.18) 

The easiest way to show this is by going to a gauge where Aq(x) is time independent. In such a gauge fl(x) = e~^ A °^ x \ 
while e~P D ° = e~^ d °e~^ A ° . But e~^ 9a — 1 due to periodicity. The equality holds in any gauge since the two operators 
e -/3D a an( ^ transform in the same way under gauge transformations. Hence, although formally exp(— f3Do) would 
be a pseudo-differential operator (Dq being a differential operator), actually it is just a multiplicative operator. 2 

The two gauge covariant constructions, [D^, ] and O(x), appear at finite temperature. The heat kernel-like expan- 
sion (expansion in powers of and X) in eq. (|2.7p is modified at finite temperature to 



(x\f{D,X)\x)=^gx(±Sl(x);T)O x (x). (2.19) 

A 

Here 0\(x) are still arbitrary local gauge covariant operators constructed with X and [D^, ]. On the other hand 
g\(±fl(x); T) are functions of the Polyakov loop and the temperature, determined by the pseudo-differential operator 
/. The ± refers to the two cases of bosonic or fermionic quantum field, respectively. 

Note that, in general, tt(x) does not commute with the local operators. We have chosen to put all the dependence 
on the Polyakov loop at the left. This can be done due to the identity 28] 

[0,g((l)]=y2-g n {il)D$0, (2.20) 
* — ' n! 

n— 1 



2 Multiplicative operators will be important in what follows. By multiplicative operators we mean zeroth order differential operators with 
respect to x. That is, operators which may contain x M but not 9 M . They can be matrices in internal space. Thus they are in one to one 
correspondence with ordinary matrix- valued functions of x. 
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where g n (£i) is just the n-th derivative of g(Q) as a function of the variable zTlog(O) and Do = [Do, }. 

For instance, the expansion in eq. (|2.19p has been computed for the heat kernel through operators of dimension 6 
infill: 

(x\e rK \x) = (ft + rtoX + ■■■), ft = £>fi) V fc2 ^. (2.21) 

At zero temperature ft = 1 and this expression reduces to that in eq. (|2.9I) . 
For the derivative expansion at finite temperature, one can write 

(x\f(D,X)\x) =^/ A (±0,X 1) ...,X n ;T)0 A , (2.22) 

A 

with the Si's at the left of all X's and 0\, and is inserted between the (i — l)-th and the i-th blocks of 0\, as 
before. (Recall that 0\ contains only operators X with derivatives. X without derivatives go into f\.) 

The functions fx in eq. (|2.22|) are well defined but have not been computed yet even for the heat kernel. This is a 
goal of this work. 

3. Countings at zero and finite temperature 

Before closing this section it is important to note that the counting of a term either by its dimension, or by its 
number of derivatives, is not as clean at finite temperature as it was at zero temperature. Indeed, to unambiguously 
classify a term by its (scale) dimension at zero temperature one can introduce a bookkeeping parameter A in the 
external fields, as XA^^Xx) and X a X(Xx) (being a the dimension of X, a = 2 in the example of eq. (|2.1|) ). In this 
way an operator O n of dimension 7 will be tagged by a factor A 7 . At finite temperature the number of X's and [D, ] 
can still be counted by a bookkeeping parameter, but the method fails for Do because a dilation in the time direction 
is not consistent with periodicity of the external fields. Of course, this is related to the presence of discrete values for 
Po and to the presence of f2(x) in addition to [Do, ]. 

At finite temperature there is no bookkeeping parameter to fix the order of a term in the dimensional expansion, 
and so the order is undefined or it may look different depending on how the term is written. To sort out this problem, 
we take the prescription of defining the counting after the term has been written with all Q(x) at the left. With this 
prescription the order can be defined without ambiguity (see appendix [A} . We take Cl(x) to be of dimension zero. As 
before X{x) has dimension a, [D^, ] has dimension one and F^ u has dimension two. So for instance, the operator 
Q(x)X(x) carries dimension a, whereas (using eq. (|2.20jl ) 

[X(x),n(x)} = f3Sl(x)[D ,X(x)}+--- (2.23) 

carries leading dimension a + 1 but is not homogeneous in this counting. As usual, we will consider the leading order 
as the order of a non homogeneous term. 

Everything is similar for the derivative expansion. In this case the zero temperature counting comes from XA^(Xx) 
and X(Xx). At finite temperature, the term is written with fi(x) at the left and then f2 and X count as order zero, 
[Dp, ] as order one and F^(x) as order two. For instance, the operator il(x)X(x) is of order zero whereas [X(x), fl(x)} 
is of order one. 

The situation for traced terms at finite temperature is more involved due to the trace cyclic property. To define 
the order of an expression in this case, the natural prescription is to consider all possible ways to write it and select 
the one with highest leading order as the true order of the expression. For instance, using the property [Do,0] = 
(from eq. (@H1), 

Tr(QX X) = Tr((QX) X) = -Tr(flXX ) = -Tr(X ftX) = Tr(-ftX X - [X 0) ti\X) 

= ~(3Tr(nX 0O X) + 0{D 3 ) = ^/3Ti(nX^ + 0(D 3 ). (2.24) 

So this is term is of second order in the derivative expansion. 

B. Symbols at zero temperature 

A convenient technique to compute the diagonal matrix elements of a pseudo-differential operator, (x|/(D, X)\x), 
is the method of symbols [13, H|| . 
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Let us discuss the zero temperature case first. The Euclidean space-time is K d 1 x R. We introduce a momentum 
basis \p), 

(x\p) = e ipx , W) = {2K) d S(p-p>) 1 |p) = e**|0), (2.25) 
and the method of symbols goes as follows 

(x\f(D,X)\x) = JJLP e - i **(x\f{D i X)\p) = J ^( x \e-^f(D,X)e^\0) 



d d p 
(2^p 



(x\f(D + ip,X)\0). (2.26) 



We have used the relations e~ tpx D ^e tpx = + ip fl and e ~ l P x X e lpx = X (because X is multiplicative, i.e., it contains 
no derivatives) and the fact that the map Y — > e~* px Ye ipx is a similarity transformation. In eq. (|2.26l) |0) is the state 
with wavefunction equal to unity, (ar|0) = 1. 

Due to the property, <9 M |0) = 0, the quantity (x\f(D + ip,X)\0) is just the symbol of the pseudo-differential operator 

f(D,X) [3a|. A very important point is that the operator J d f(D + ip,X) contains only in the form [D^, ]. 

As a consequence, this operator is automatically gauge covariant and also multiplicative with respect to x. As said, 
a multiplicative operator is equivalent to a function of x. Specifically, (x|/(x)|0) = f(x)(x\0) = f(x). So (a;| |0) can 
be left implicit in eq. (|2.26p . and one can write just 

(x\f(D,X)\x) = J jULf(D + ip,X). (2.27) 

The variable p M represents the momentum carried by the quantum field <j> running in the loop. 

To obtain a covariant derivative expansion, one simply expands the right hand side of eq. (|2.27p in powers of D^. 
Due to gauge invariance, it is guaranteed that if all are brought (e.g.) to the right using D^Y — Y] + YD^, 
at the end all terms with not in the form [D^, ] must vanish after momentum integration. So gauge invariance of 
the final result will hold but it is not manifest without momentum integration. 



C. Covariant symbols at zero temperature 

The matrix element (x\f(D, X)\x) is a gauge covariant quantity, and its covariant derivative expansion can be 
obtained by expansion in powers of in eq. (|2.27[) . However, gauge covariance of the right hand side holds only 
after momentum integration: the symbol itself is not covariant. Plctnev and Banin devised a method to transform 
the symbol into a covariant one [HI, [42| . This is as follows 

(x\f(D,X)\x)= J -0^f(D,X). (2.28) 

with the covariant symbol 

f(D, X) = e ld " D e- ipx f(D, X)e lpx e~ tdV D , d p = — : d p D = Dd p = D,d p (2.29) 

That is, a further similarity transformation is applied which changes nothing: the new factor e~ ldFD is equivalent to 
1 since no p^ lies at its right, and on the other hand the new factor e ld D is also equivalent to 1 by integration by 
parts. Being a similarity transformation it can be applied to each block in /, i.e., — > and X — > X. 

D„ = e^ D e- ip *D^ p *e- idPD =e id ° D {D IM + ip li )e- idPD , 

X = e i9 p D e ^px Xel px e - t d"D _ e id*D Xf ,-id>>D ^ ^ 

These new operators are directly gauge covariant and multiplicative (with respect to x) without momentum inte- 
gration. Using a derivative expansion, they read: 

oo 

n.—l ^ '' 

oo 1 

1 = £^ <BJf «i-«nflS 1 -SS n - (2-31) 

n=0 
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As can be seen, the covariant symbol is closely related to the Fock-Schwinger gauge approach. The map Y —> Y is an 
algebra homomorphism that applies pseudo-differential operators into operators which are covariant and multiplicative 
(with respect to x). They are derivative operators with respect to p^. Let us stress that, in applications of eq. f|2.28j) . 
a constant function equal to 1 is understood at the right, so that 1 = 0. 3 

D. Symbols at finite temperature 

Let us now turn to the finite temperature case. For ordinary symbols one can proceed as before by introducing 
a momentum space basis \p) = \po,p), where the zeroth component takes values on the Matsubara frequencies: 
Po = 2irnT in the bosonic case, po = (2n + l)wT in the fermionic case, with n 6 Z. Thus 

(x\p)=e i ** ) {p\p')=^ Po , p ^) d - l 8[p-p'), |p)=e**|0). (2.32) 

The method of symbols works as before with the result 

(x\f(D,X)\x) = tW ^-^ T (x\f(D + l p,X)\0). (2.33) 

Po \ 7T> 

Let us remark that |0) is the state (acjO) = 1, regardless of whether the quantum field in the loop is bosonic or 
fermionic. The statistics of the quantum field is contained in the Matsubara frequencies po- Once again the operator 
TJ2 Po J ( 2 °7<i-i f(D + ip,X) is actually multiplicative and (x\ |0) can be omitted 

(x\f(D,X)\x) = tW ^^f(D + ip,X). (2.34) 

Po 

Also (x\f(D, X)\x) is still gauge covariant. 

In previous works we have discussed the effect of the finite temperature, i.e., the replacement of an integral over 
Po on R to a sum of po over Matsubara frequencies. As in the zero temperature case, after integration over p, the 
operator D appears only in the form [D, ]. The reason is that obviously if one replaces D by D + ia, a being 
a constant c- number, the replacement has no effect owing to the integration over p on R''" 1 . However, the same 
argument fails for Do (the zeroth component of the gauge covariant derivative) since po is a discrete variable at finite 
temperature. Still, due to the sum over the Matsubara frequencies, the expression must be periodic in the variable 
Dq with period 2niT. This not only permits a dependence on [Do, } but also on e^ 130 " — J7, i.e., on the Polyakov 
loop. 

Let us discuss how to use the ordinary symbols to obtain the diagonal matrix elements at finite temperature [33j | . 
The main issue is the gauge invariance. In the method of symbols, eq. (|2.34j) . gauge invariance of (x\f(D,X)\x) 
is manifest only after the integral on p and the sum on po are carried out. In f(D + ip,X), D can be dealt with 
as in the zero temperature case to yield \ D, } after integration on p. This produces an expression of the type 
fi(D +ip ,[D, ],X). As described in [HI, |33|, a method suitable to deal with D to obtain a derivative expansion, 
is to move Do to the left (using the identity YDq — DqY — [Do, Y]). In this way one ends up with expressions of the 
type feiDa + Wo] [Do, ], [D, ],X) where Dq + ipo is only at the left, rather than all over the expression. Summing 
now over the Matsubara frequencies produces a dependence on e^^ D,) = i7 and finally a covariant expression of the 
type h(Vt, [Dp, },X) with all Q,(x) at the left. This is the form in eq. (|2~19| or in eq. (f2T22]> . 

E. Covariant symbols at finite temperature 

The method just described at the end of the previous section is rather cumbersome, so a method of covariant 
symbols at finite temperature would be advisable, namely, a method providing manifestly multiplicative and gauge 
invariant terms. The problem is that the method of Pletnev and Banin does not directly apply at finite temperature, 
since po is a discrete variable and <9q is not defined. Presently we show how to extend the method to the finite 
temperature case. 



3 Actually, the quantity f (x,p) = f(D, X)l is what enters in the computation of (x\f(D, X)\x). This is an ordinary function of x and p 
and so closer to the ordinary symbols, except that it is covariant. 
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One idea is to change the sum over Matsubara frequencies by appropriate integrals on the complex plane [22J. In 
this way the derivative with respect to po is defined. This method works but we can obtain the final result in a simpler 
manner. 

Let u) n be the Matsubara frequencies, bosonic (u;„ = 2mtT) or fermionic (w n = (2n + l)irT). Then let 

h M (Po) = ]T 2nTS( Po - u n ). (2.35) 

n 

(There is a bosonic version and a fermionic version of this function.) 
Using the function Km we- can write eq. (I2.33|) as 



d d p 
{2tt)° 



(x\f(D,X)\x) = I 7 *h M (jp )(x\f(D + ip,X)\0). (2.36) 



Now we can proceed to make a further similarity transformation, as at zero temperature (and valid by the same 
reasons) 

(x\f(D,X)\x) = J ^( x \e* dPD h M (po)f(D + l p,X)e-> dPD \0) 

' l '' l> (x\e tdPD h M ( Pa )e- ldFD f(D,X)\0). (2.37) 



(2w) d 

This can be simplified by working out the /im(Po) term: 

e WPD p e- w " D = Po + iD - ±iF 0i df + ^flgflf + ■■■ (2.38) 

hence 

e WPD h M ( Po )e- tdPD = h M (po + iD ) + 0(9?). (2.39) 

The point is that, due to the integration on p, all df at the left (no p lies at the left of the df) can be set to zero, 
and so 



d d p 
(2nj 



(x\f(D,X)\x) = I 7 ^ 3 {x\h M (po + iD )f(D,X)\0). (2.40) 



The expression eq. (|2.40l) is of great interest. D and X are the same covariant symbols as at zero temperature, 
and so they are Lorentz covariant (if the original pseudo-differential operator / is). They are also multiplicative 
with respect to x-space and manifestly gauge covariant. On the other hand the Dq dependence at the left is also 
multiplicative: under the shift Dq — > Dq + 2mnT the expression is unchanged due to periodicity of Km (even without 
integral over po). Therefore, the dependence is really on the periodic variable e~^ Da = fl. That is, one can also write 4 

(x\f(D,X)\x) = J ^h M (p -iT\o g n)f(D,X). (2.41) 

This expression is already of the form required, gauge covariant and with at the left, suitable to take the expansions 
in eq. ([2~T9l> or eq. (|2~22l . 

For convenience let us introduce the auxiliary multiplicative operator (a matrix in internal space) 

Q(x) =iT log fl{x). (2.42) 

This is many- valued but in practice it appears in periodic functions so that the result is always a one-valued function 
of 17. Q is Hermitian, up to many-valuation, fl being unitary. Eq. (|2.4ip takes the form 

(x\f(D,X)\x}= J *Ph M (po-Q)f{D,X). (2.43) 



4 Once again, in eq. 112.411 1. a constant function equal to 1 is implicit at the right, so that 1 = 0. 
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It is possible to define also the quantity Qq as the operator Q placed at the left of all other operators, that is, 
labeled to indicate "at position zero" . There can be no confusion with our previous convention of a label indicating 
a temporal covariant derivative since [Do,Q] — due to [Do,0] = 0. The point is that Qq is a c-number. it can be 
put in any order in an expression with the same result. Hence we can shift the variable po by an amount Qq. This 
allows us to write 

( x \f(D,x)\ x) =TJ2J |^/(Ax)|^ +Qo 

- T E/ (Si/^o+^A,!). (2.44) 

po 

(In the last equality we have used that the variable po does not appear in Dq — ipa, Di or X.) 

In eq. (|2.43l) one can carry out the momentum derivatives d p implied by and X. The derivatives df can be 
taken to the right or to the left, by parts. The temporal derivative 8q can only be taken to the right, if the form 
of /im(.Po — -Do) is to be preserved. Taking all of the d p to the right has the virtue of leaving an ordinary function 
f'(x,p) which is temperature independent, and manifestly Lorentz and gauge covariant: 

(x\f(D,X)\x) = J APh M (po-Q)f'(x, P ). (2.45) 

The eqs. (|2.43[) or (|2.45[) solve the problem of using gauge covariant symbols at finite temperature. In addition the 
breaking of Lorentz covariance is minimal. The zero temperature limit is recovered by setting Hm to unity. 



F. Polyakov loop and real time thermal field theory 



Mathematically, the imaginary time formalism is the simplest approach to quantum field theory at finite tempera- 
ture. The real time approach (in its various versions) is more involved but better suited for time dependent observables 

In that approach the frequency is a continuous variable, rather than discrete. In the expressions derived for the 
covariant symbols at finite temperature the sum over Matsubara frequencies can be traded by integrals using well 
known relations (eqs. (2.3.22-24) of 22]). Starting from eq. (|2.44p . we find 

(x\mX)\x) = J ^L(||!l + ^n(- PQ )-J c ^-n( P0 ))f(D + i Q ,D l ,X). (2.46) 

Here 

»<*>) - T^—r (2 - 47) 

(± for bosons or fermions, respectively). The first frequency integral is along the po real axis, whereas the contours 
C± enclose only the singularities of / as a function of po, in the half planes Impo > and Impo < 0, respectively. 
Undoing the shift po — > po + Qo, gives 

(x\f(D,X)\x) = [ f~*P ([**+[ *° -J-, / *° i )f(D,X). (2.48) 

V UV ' 71 1 J (27r) d - 1 \J 2vr J c+ 2tt ifi-^-" 9 ™ - 1 J c _ 2ir ±Q,e^P° - l) J V ' ' V ' 

(C± are as before for the new /, since Qo is real.) 

This expression is not yet in the form of the real time formalism but it is closer to it. Upon Wick rotation, factors 
of the type n(pg) should appear there in the propagators through the thermal occupation numbers, while the integral 
over the real axis should come from the zero temperature part of the thermal propagators [22| . 

The connection with the real time formalism is, of course, of great interest, and worth studying. We do not pursue 
this subject any further in the present work, but in view of eq. (|2.48p one can conjecture that the fields in the form 
f(D,X), including time covariant derivatives, will follow the pattern of ordinary local external fields as treated in 
the real time formalism. On the other hand, the occupation number will pick up a Polyakov loop following the 
prescription e@ Po — > fle^ Po . This automatically produces the correct coupling of the chemical potential, e^, by means 
of the prescription Aq(x) — s- Aq(x) — fi (since — Aq essentially represents a local and possibly non Abelian chemical 
potential) . 
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The meaning of Cl(x) in the real time context needs to be elucidated. If the configuration of the external fields is 
stationary, essentially the imaginary time formula gives already the result. In this case fi(jc) = e ~ pA °^ with Aq(x) 
Hermitian. Aq is the same variable in the Euclidean and Minkowski versions, however, in the Minkowski case this 
variable is taken along the real axis of its complex plane, while in the Euclidean version it is preferable to work with 
the variable extended to the imaginary axis. The imaginary time integration in the definition of is not rotated to 
real time because it comes from the factor e~^ H in the par tition function and so from an evolution in imaginary time 
from to to to — i/3 (Kubo-Martin-Schwinger condition [43|). Therefore takes the same form any finite temperature 
formulation. 

Let us consider now the more general case of non stationary configurations. In the closed-path approach [44|,|45| one 
starts with a thermal mixed state at time t — t . This implies that the system is stationary for t < to- Measurements 
are taken at later times, where also time dependent sources may act. By assumption of thermal equilibrium, either 
the Hamiltonian is stationary for t < to, or it is so upon taken a suitable gauge transformation. Then, for t < to, 
Aq(x) is well defined modulo stationary gauge transformations and this defines fl(x) — e~^ A °^ which also transforms 
covariantly. Further gauge transformations, for instance carrying Aq to zero, exist but they are not stationary and 
so they would introduce a time dependence in the other components of the gauge connection (and possible on other 
gauge covariant fields). Therefore such transformations are not allowed for t < to and f2 is well defined. Because the 
external fields configuration is not required to be stationary for t > to one can choose a continuous gauge in which 
Aq(x) takes the same value at all times. (Because any Ao(x,t) can be brought to zero by means of a suitable gauge 
transformation, any configuration Ao(x,t) can be transformed into any other.) This shows that Q(x) is also present 
in the real time approach and this is the quantity that will appear with e^ Po in the propagators. 

III. HEAT KERNEL AT FINITE TEMPERATURE 
A. Diagonal coefficients 

1. Expansions of the heat kernel 

Let K = D 2 + X be the Klein-Gordon operator as in eq. (|2.ip . The heat kernel is the solution of the associated 
heat equation d T G{r) = KG{t), G(0) = 1, r > 0, with solution G(r) — exp(rif). 5 From the heat kernel one can 
recover the propagator, A" -1 , and the effective action, Trlog A. 

The diagonal matrix elements of the heat kernel (at zero or at finite temperature) can be expanded classifying the 
terms by their mass dimension: 

{Ae rK \x) = j^ w Y t r-a n {x-,r), (3.1) 

Each a n has dimension In and depends on the temperature. The expansion is asymptotic. At zero temperature this 
is equivalent to an expansion in powers of r, and this is just the standard heat-kernel expansion. In general the a n 
depend also on r and T. The order of the term is defined by the mass dimension carried by the external fields. Hence 
by dimensional counting, the coefficient can only depend on the combination tT 2 . A remarkable property of the heat 
kernel coefficients is that they do not depend explicitly on the space-time dimension. This property is preserved at 
finite temperature. 

At zero temperature the index n takes non negative integer values. 6 However, at finite temperature n can also 
take (positive) half-integer values. This follows from breaking of Lorentz invariance down to rotational invariance; at 
finite temperature an odd number of time derivatives is not forbidden. The expansion at finite temperature has been 



5 For the Klein-Gordon operator the parameter r has dimensions of inverse mass squared, nevertheless it is called the Fock-Schwinger 
proper time [46l| since in the heat kernel equation it plays the role a time with corresponding Hamiltonian iK acting in the Hilbert space 
spanned by \x). 

6 There are half integer orders in the presence of boundaries. We only consider boundaryless manifolds throughout. 
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computed in [22],[28[ through dimension 6. So for instance, 7 

ao = 6, 
a-i/2 = 0, 
a\ = £oX, 
1 
2 



« 3 /2 = -6(*o+£«)- (3.2) 



The electric field, Ei(x), is defined as Foi(x), hence En — —Fuq. On the other hand, the 6 are dimensionless 
functions of the Polyakov loop defined as sums over the (bosonic or fermionic) Matsubara frequencies: 



6 = (4^r) 1 / 2 (-*)"2-"/ 2 r^i7 n (V27(po + 0))e- r(po+C)r ' 

Po 



= 2-"/ 2 ^i/„(/fc/v / 27T2)e- fe2 /( 4rT2 )(±0) fc , n = 0,1,2,... (3.3) 

feez 

Q was introduced in eq. (I2.42p . H n refers to the n-th Hermite polynomial (with normalization H\{x) — 2x). The ± 
refers to bosonic or fermionic case, respectively. The two forms of 6 in eq. (|3.3j) are related by Poisson summation 
formula. The 6 are one- valued functions of f2 and of tT 2 . They are real (Hermitian) for even n and imaginary 
(antiHermitian) for odd n. In addition, they are even or odd under il —> SI -1 for even or odd n, respectively. In the 
zero temperature limit 

e x=o = 2 -»/2fl n (o) 1 (3.4) 

so odd orders vanish in this limit. 

It will be also convenient to define the following auxiliary combinations 

6 = 6, 6 = 6+ Co, 6 = 6 + 36, 6 = 6 + 66 + 36- (3.5) 

They vanish at zero temperature. (However, the 6 do not vanish at finite temperature for Q = 1 for even orders.) 
The derivative expansion of the heat kernel (at zero or finite temperature) takes the form 

^ TK \x) = J^ynY,- nA ^T), (3.6) 

where the coefficient A n contains 2n derivatives, as well as the Polyakov loop (placed at the left) and any number of 
X. By dimensional counting, besides the derivatives, A n (x; r) depends on tX and tT 2 and Q. This is an asymptotic 
expansion. Once again, at zero temperature the index n takes only nonnegative integer values, whereas at finite 
temperature half-integer values are allowed. The derivative expansion coefficients A n are also independent of the 
space-time dimension, at zero or finite temperature. 

The expansion at zero temperature has been considered in [35j to four derivatives (and six derivatives for the traced 
coefficients). For instance, 

A) = h, 

Ax = Tl 2a X^ + 2T 2 I 2X2 Xl. (3.7) 

The coefficients Ii, 7 2 ,2 and 1% 1,2 are functions of the labeled operators X\ in the first case, X\, X 2 in the second, 
and X\, X 2 and X 3 in third case. In general, these coefficients are defined as follows [35j 

Ir u r 2 ,...,r n = J ^e?N? N? • ■ ■ N?, n = 0,1,2,... (3.8) 

where 

N t = {z-rXi)- 1 , (3.9) 



7 Regarding conventions, let us note that what is called here K and X corresponds to — K and — M in [27l 1 28t] ■ The functions £ n are 
similar to the ip n in [27l . |28| except that they involve the Hermite polynomials. 
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and T is a positively oriented simple closed path enclosing all the AV 8 Explicitly 

I - T i-E?=i n V - - - (3 10) 

The functions I ri ,r 2 ,...,r n are analytical on the Xi even at coincident points (as follows from eq. (|3.8j> . the singularities 
at Xi = Xj are removable) and satisfy recurrence relations. Instances at lower orders are 

e rXl 

Ir = t -T7, r = 0,1,2,... 

(r- 1)! 



1 e rXl + e rX2 2 e rXl - e rX2 

7* (x l - x 2 y ~ 73 (X l - X 2 f 



J 2.2 - ^Y2~Z3— V^W' ( 3 - 11 ) 



2. Derivative expansion at finite temperature 

The coefficients A n at finite temperature are not yet known. They can be computed from scratch by using the tools 
previously described. To this end we use an integral representation of the heat kernel 

where the path T is positively oriented and encloses the eigenvalues of K (the concrete realization of this requirement 
will be clear below). 

Applying the method developed in section III El for covariant symbols at finite temperature, and in particular 
eq. (|2.43p . we can write 



Using the explicit expressions of the covariant symbols of and X in eq. (|2.3ip . it is simple to carry out an 
expansion with terms classified by the number of covariant derivatives they have (regardless of the number of X or 
Q). Specifically, 9 removing the zeroth order contributions in L> M and X, 

D'^ = D^-i P ^ = 0(D 2 ), X' = X-X = 0(D), (3.14) 

we can write 

(z-Di-X)- 1 = (N-'-iip^D'^-D'Z-X')- 1 

oo 

= Y l N((i{p„D' ll }+D*+X t )N) n , (3.15) 

n=0 

where we have introduced the quantity 

N = (z+p 2 -X)- 1 . (3.16) 

Let us spell out the details for A\i<z (i.e., one derivative). Picking up the terms with precisely one derivative in 
eq. (|3.13[) gives (using eq. p.!5|) and eq. (|2.31jl ) 

(x\e TK \x) 1/2 = JJ^J ^^(po - Q)e r 'NiX M d^ t N. (3.17) 
Further, applying the identity 

(d?N) = -2 Vii N\ (3.18) 



8 This r is not to be confused with the effective action functional introduced in eq. I|2.2t . 

9 Alternatively one can use the formulas of Appendix \C\ for the covariant symbols of K and [z — K)~ 1 . 
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yields 



(x\e rK \x) 1/2 = ]J^J 4^d h M(Po ~ Q)e TZ (-2i)p li NX li N 2 . (3.19) 
Next, let us apply the shift z — > z ~ p 2 , so that 

(x\e TK \x) 1/2 = J ^Lh M ( P0 - Q)e-^\-2i)p^ J^e^NX^N 2 , (3.20) 

where 

N = (z-X)~ x . (3.21) 
Now the z and p integrals are independent. For the z integral, the definition of I ru ...,r n in eq. Q3.8[) applies: 

^-.e TZ NX^N 2 =T 2 h, 2 X^. (3.22) 
Zttz 

For the p integral, the definition of £„ in eq. (I3.3[) applies: 

" lp jh M (p Q -Q)e-^ 2 (-2 l )p^ = (-21)5^ - [ ^h M (p - QJe^po 



(2ir) d mKVU ^' y ^ ~ v ' Aiu (4 7 rr)(' i -i)/ 2 J 2ir 

1 

(47TT) d / 



"MO ,, \d/2 T ~ 1/2 fr- ( 3 - 23 ) 



Therefore, 



<*|e T *l*>i/2 - J^^r^^Xo. (3.24) 

or according to eq. (j3.6j> . 

Ai/a = rCi/i,2^o- (3.25) 
In what follows we use units r = 1. t can be easily restored by dimensional considerations. 

Using the method just described and the formulas in Appendix [B] for the momentum integrals, we find to three 
derivatives 

A = iifo, 

Ai/2 = h,2^lXo, 

A\ = 72,2 £o X^ + 2/2,1,2 Co X^X^ + 71,3 £2 ^00 + (2Ii,l,3 + -^1,2,2) £2 XqXo, 

1 - 2 - 

A3/2 — h^^iFofiXfj + -Ii^iFftOfj, + -/2,3d (-^0^ + X^o^ + X^o) 

+ (2/1,1,3 — 6/14,4 + /i,2,2 — 2/1,2,3) S.iXqX^^ + 2/2,1,3^1 (X^Xq^ + X^X^q) 

+ (2/2,1,3 + /2,2,2) £l (-X'o/j-X'/i + XftoXf, + X^^Xq) 

+ (4/i,2,i,3 + 2/1,2,2,2 + 4/1,3,1,2 + 4/2,1,1,3 + 2/2,1,2,2 + 2/2,2,1,2) S.iXoX^Xn 
+ (4/2,1,1,3 + 2/2,1,2,2) ^{X^XoXfj, + X^X^Xo) + /1, 4^3^000 

+ (3/l,l.4 + /l,2,3j ^3^0-^00 

+ (3/l,1.4 + 2/i,2,3 + h,3,2) ^3^00-^0 

+ (6/l,l,l,4 + 4/i,i,2,3 + 2/i,i,3,2 + 2/i,2,l,3 + /l,2,2,2j ^3^0^0X0- (3.26) 

Important remark: For notational convenience we have written £0 or £n at the right of the I ri ,r 2 ...., ou t actually these 
operators are at the left of the expression. So Aq = £0/1, -A 1/2 = £1 h^Xo, A\ — ^0/2,2-^^ + • • • % etc. 

We have also computed the term with four derivatives A 2 , but this term is too long to be quoted here (about 90 
terms). The four derivative term is given below for the traced heat kernel coefficients. 

The heat kernel (and in fact (x\f(K)\x) for any f(z)) is symmetric under left-right transposition of operators (or 
Hermitian if K is Hermitian and f{z) is real). At zero temperature (putting — > and £0 — > 1) the symmetry 



15 



is manifest. For instance, the term + 2/2,1, 2X fl X IJi is symmetric. The symmetry is not manifest at finite 

temperature because it is hidden after having chosen to put the Polyakov loop to the left. In addition, transposition 
and subsequent move of the Polyakov loop to the left in a term A n produces new terms of higher order. For instance, 
to first order in the derivative expansion, 

(£0/1 + £Vi,2*o) T - + £1/1,2*0) = [A, Co] - Zi(h,2 + h,i)X + 0(D 2 ). (3.27) 

From cq. (|2.20[) . [/i,£o] = [-Do, h]- Use of — = [-D ,/i] = h,iX , and / li2 + /2.1 = h,i, shows that the 

symmetry holds to the order considered. 

The /ri,...,r n are not linearly independent, so although the coefficients in A n are well defined functions of the labeled 
operators, X,-, their expression in terms of the I ru ... jTn IS n °t unique. 

The heat kernel does not depend on the prescription adopted regarding the position of the Polyakov loop (our 
choice throughout has been to put it at the left) but the value of each A n will be different for different prescriptions. 

The Polyakov loop comes out automatically in the expressions, and, as noted in the Introduction, its presence is 
required to accommodate the chemical potential. Nevertheless, it is also a nuisance and so the possibility suggests itself 
to dispose of the Polyakov loop dependence just by setting O = 1 in the formulas by hand. We call this the quenched 
approximation. If this is done, the become ordinary functions of the temperature (rather than operators) for even 
n and zero for odd n. Unfortunately, the result of quenching will depend on the prescription adopted (regarding 
the position of the Polyakov loop) and in particular the left-right symmetry can be lost. In fact, the expressions are 
fully consistent only when the full Polyakov loop dependence is retained. For the traced heat kernel, and so for the 
effective action, setting 51 to unity by hand is also dangerous. Due to the cyclic property, the same expression can 
be written in several equivalent but different ways. Consequently, the result obtained by setting — 1 by hand will 
yield different results in each case. This point is further discussed in Section HV C 41 

B. Traced heat kernel coefficients 

It is also of interest to compute the trace of the heat kernel and this produces shortest expressions. Specifically 
(remember that we have set t = 1) 



Tr(e^) = J d d xtv(x\e K \x) = ]T J d d xtrB n (x) 



(3.28) 



The choice B n — A n is of course correct, but some simplification in the form of the coefficients B n can be achieved 
by using integration by parts and the cyclic property of the trace. When using this freedom, the functions should 
be moved to the left by using the identity in eq. (|2.20[) . 

Note that the A n can be recovered from the B n using the identity 

This equality holds separately at each order in the derivative expansion. 

For convenience, we separate in B n terms with a contribution at zero temperature from those which vanish in that 
limit, 



B n = B^+BlP. (3.30) 



The Bn^ vanish for half-integer n, and are of the form £o(-S n | T _ ), while B^\ TQ = 0. 10 
The results are as follows: 



10 Note that the definition B^ = £o(£?n| T _ ) would be ambiguous, since B n \ T _ can be written in different ways which are equivalent 
inside the trace, but not in £,a{B n \ T _ Q ). 
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(0) 



D. 



(o) 



hto, 

-■^h^^oX^X^, 

+^h,i,3,o^oXf 1 X li X I/v + -ha^afiioX^XyX^X^ 
+(4/3,1,3,1,0 — hfi^fifl^oX^X^X^Xv. 



(3.31) 
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B (T) - 

b[ T ^ = -^h,2,i&XoX , 



B 3/2 = (-^^1,2,0-^-^2,1,0)6^^ 



1 , 1 



g -^1,2,2 - g ^1,3,1 ) ( Cl ^0,, A ,, - i A ..... A A : A ; . /; A 1) - : -i„; A (ID ^ i) 



( 2^1,1,2,2 



3/1,1,3,1^ (^lXoX^Xfj + ^X^XqX^ + ^X^Xf+Xo — -£3^o^o^o^ 



+ ^3-77/3,2,0 - 2^ 3 ' 30 ~ 2^ 4;2,!l ) 7 

gg ^1,3,0 — 3/2,2,0 — 3g^ 3 > ] ' ) ^--A ( i/i ' (),/ 

'7 1 1 

^^3,2,0 _ TJ^ 3 ' 3 ' _ 2^ 4 " : '" ! ' V ''' ,A "" 

1 , 1 



+ ( - 2^ 3 - 3 ' _ 2^ 4,2,0 / (^ 2 ^°^^° + ?2^mo^om + £2-^0^,0 - 2^ 4 ^ 00 ^ 00 
7 11 \ - 

+ ( ^2,1,3,0 — ^2, 1,4,0 + 3/2,3,2,0 + JT^^, 1,2,0 — 2/3,1,3,0 — I4, 1,2,0 j ^XqXqX^^ 



+ 



11 1 13 1 11 17 

3^/1,1,3,0 - y^i.2,2,0 - 3^/1,3,1,0 - 3/2,1,2,0 - 3^/2,2,1,0 - 3^/3,1 

11 5 11 1 1 17 

^/l, 1,3,0 + -^7/1,2,2,0 — ^/l,3,l,0 — o/2,l,2,0 — 77/2,2,1.0 — -077/3,1,! 

36 36 3d 3 9 36 

7 17 \ - 

gh, 1,3,0 - h, 1,4,0 + 3/2,3,2,0 + g/3,1,2,0 - 2/3,1,3,0 - /4,1,2,0 J (^2^0^^0m + &X ^XqXqh) 

1 17 \ - 

Y^/2, 1,3,0 — h, 1,4,0 + 3/2,3,2,0 + y^/s, 1,2,0 — 2/3,1,3,0 — h, 1,2,0 1 £2 "^u^/^oo 

+ ^-/2, 1,4,0 + 3/2,3,2,0 - 2/3,1,3,0 - I4, 1,2,0^ ^-Xb-Xfi-^O + £2-Xf*-Xb-Xft0 - 2^ 4 ^0^0^00^ 

+ f jtj^ 2 ' 2 ' 2 ' 2 ' g/3,1,2,1,0 — 2/3,1,3,1,0 + 3/3,2,1,2,0 — 2/4,1,2,1,0 J ^XqXqX^X^ 



5 7 2 \ - 

Y^/2,2,2,2,0 + g/3,1,2,1,0 — 2/3,1,3,1,0 + 3/3,2,1,2,0 - 2/4,1,2,1,0 J i&XoXpXo Xp + £2X0X^X^X0 

+(,2X^XqXqX^ + £2X^X0X^X0 + £2X^X^X0X0) 

+ ^—7^/2,2,2,2,0 + /3, 1,3, 1,0 — 3/3,2,1,2,0 + I4, 1,2, 1,0^ £4-^0-^0-^0-^0- (3.32) 
Once again we note that the £o and £n are actually at the left of the / ri ,..., r „- 

Further rearrangement of the expressions is possible to bring them to a more systematic form. For instance 
reordering of covariant derivatives is possible using the Bianchi identity Y^ v = Y vfl + [F^ V ,Y] 7 as well as cyclic 
permutations or integration by parts. However, such extra work does not seem to yield a simpler expression. These 
expressions for B n have not been obtained directly from A n but from Tr log(z — K) in Chan's form, to be introduced 
below. 
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IV. CHAN'S FORM OF THE EFFECTIVE ACTION 



Up to now we have considered Euclidean space-times with the topolo gies R d or x S 1 appropriate to study field 
theories at zero or finite temperature. The latter case leads to the Matsubara frequencies and to the weight function 
1im{Pq) introduced in eq. (I2.35|) . At zero temperature the weight function is just equal to unity. 

As it turns out, the formalism can be carried out equally well without assuming any particular properties of the 
weight function h(p) in the momentum integration. h(p) can even depend on all components p^. For the purpose of 
deriving general expressions no simplification is obtained by imposing constraints on h(p), therefore, from now on we 
will assume a completely general weight function h{p). We call /i-space the setting leading to such a weight h{p) in 
the momentum integrals. In the next subsection we show that this approach does not lead to inconsistencies. 



A. /i-spaces 



We devote this subsection to study the consistency of the approach with generic h(p), specifically regarding gauge 
invariance and cyclic property. 

Generalizing the method of symbols, we define 

(x\f(D,X)\x) h = J ^h(p)f(D + ip,X), (4.1) 

Tv h f(D,X) = Jd d xtv(x\f(D,X)\x) h = J^fh{p)tvf(D + ip,X). (4.2) 

h(p) is a c-number function, therefore the cyclic property works as always: Tr^(/i/2) = Tr^^/i). 11 As a conse- 
quence the following property holds 



STr h (e K ) 



5X{x) 

To extend the method of covariant symbols for generic h(p) we define 



= (x\e K \x) h . (4.3) 



h(p + iD) = ]T -ffid^ ■ ■ ■ dl n h{p))D, n D^ ■ ■ ■ D^. (4.4) 

71=0 

Then 

(x\f(D,X)\x) h = J ^ h{p )e-^e iDdP f(D + i P ,X)e- iD9P 
d d p 



(27T) 



h(p + iD)f(D,X). (4.5) 



Let us consider now the issue of gauge invariance of (x\f{D,X)\x)h- To study this, it is convenient to write the 
r.h.s of eq. (|4.1[) more explicitly as 



d d p 
(27r)° 



(x\f(D,X)\x) h = I 7 -JLh(p){x\f(D + ip,X)\0). (4.6) 



Now, any operator O constructed with and X(x) necessarily transforms gauge covariantly, i.e., as U~ 1 OU. 
Gauge covariance can be lost by taking matrix elements with the state |0), which is not covariant: in general (x|0|O) 
does not transforms into U~ 1 {x){x\0\Q)U(x). However, the correct transformation is guaranteed provided O is a 
multiplicative operator because in this case (x|0|O) = O(x)(x\0) = 0(x) —> U^ 1 (x)0(x)U(x). Therefore, gauge 
covariance of (x\f(D, X)\x)h is ensured provided the operator 

/' = / ^h(jp)f{D + ip,X) (4.7) 



1 When h(p) = 1, a good convergence of f\ 2(D + ip, X) for large p M is assumed. Here we assume that this convergence is not spoiled by 

Hp). 
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is multiplicative (matrix elements (x| |0) have not be taken here). The same requirement holds for the operator 
h(p + iD) in eq. (|4.5[) . namely, it must be multiplicative. (The covariant symbol f(D,X) is already gauge covariant 
and multiplicative.) 

For an operator O to be multiplicative amounts to commute with c-number functions of x. This requirement can 
be recast in the form (the fc M are constant c-numbers) 

e -ikx Qe ikx = Q ^ 

Due to the property e~ %kx D ^ kx = + ikp, we can see that /' or h(p + iD) will commute with e zkx if 

h(p-k) = h(p). (4.9) 

If this condition is imposed for all k, the function h(j>) must be a constant. This corresponds to the zero temperature 
case. In this case, the quantum fields belong to the vector space Va of arbitrary functions of x in R d (we disregard 
internal degrees of freedom here). At finite temperature, the quantum fields arc required to be periodic or antipcriodic, 
and the external fields periodic. This implies that one is working now in a subspace V of Vd, (namely, that of periodic 
or antiperiodic functions). The operators (external fields) acting on that space can carry only momenta of the type 
k = (fe,cj„) in order to leave V invariant. Therefore, one needs to consider only this set of momenta when checking 
the relation h(p — k) = h(p) for h(p) = hu{po) (and the relation is of course fulfilled by 1im(po)-) At the same time, 
the restriction in k is directly related with the compactification M. d — > x S 1 . 

Let us generalize these ideas for other h{p). There should be a vector space V of space-time functions for the 
quantum fields, a set A of allowed operators leaving V invariant, and a set K of allowed momenta. The operators in 
A are those having only momenta A; in if in their decomposition in Fourier modes. Because combinations of operators 
in A should also stay in A, we must demand that if k±, ki £ K, k\ ± k2 € K (i.e., the set K is closed under linear 
combination with integer coefficients, in particular € K). On the other hand, V is composed of those functions 
with Fourier modes of the type q + k, for some fixed q and k £ K. Ideally, such K would come from some suitable 
compactification of R d . Finally, there will be gauge invariance provided h(p — k) = h(p) for all p and all k in K . 

In practice the only obvious setting carrying out the above program is for space-times of the type 1" x S 1 x • ■ ■ x S 1 , 
< n < d. This corresponds to modes k which are an integer linear combination of d — n fixed linearly independent 
vectors plus an arbitrary vector in the n supplementary directions. In this case h(p + iD) is a function of the d — n 
"Polyakov loops" in the d — n compactified directions. 

At present, it is not clear whether there exist other useful realizations of /i-spaces. In any case, the formalism can 
be developed without special assumptions on h(p). In what follows we simply assume that the quantum fields lie in 
the appropriate space V (the /i-space) and the allowed external fields, as well as the allowed gauge transformations, 
leave V invariant. 



B. X-form and iV-form of the expressions 



1. Diagonal matrix elements of the propagator 



Let the propagator be 



G(z) = (4.10) 



As is well-known one can obtain generic functions of K from the propagator, 

f dz 

f{K) = j v i£-f{z)G{z), (4.11) 

where T encloses counterclockwise the spectrum of K. (f(z) is assumed to have the required good properties.) 

The diagonal matrix elements of the propagator in the /i-space can be computed using the method of symbols or 
covariant symbols and the derivative expansion, as already explained for the heat kernel. To second order one finds: 

(x\G(z)\x) h = J -jL*h{p + iD)(N -2ip^NX^N 2 - Ap, lPu NX^N 3 + NX^N 2 

~Sp^ Pu NX^NX y N 3 - Ap^p u NX^N 2 X v N 2 + 2NX^NX, L N 2 + 0{D 3 )) . (4.12) 
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Here 

TV = (z+p 2 -X)- 1 . (4.13) 
The expression through third order is given in Appendix [C] using labeled operators. 

We refer to the form in eq. (|4.12l) as the X-form of the expression because the X appear with derivatives and the N 
carry no derivative. By means of the relation X^ — N~ 1 N fl N~ 1 , and derivatives of it, one can eliminate completely 
the X and write the same expression using only N and covariant derivatives of it. For a generic initial expression, 
negative powers of N will be present after elimination of X . When this is not the case we say that the expression 
admits an N-form. As it turns out, the covariant symbol of the propagator admits an iV-form (see Appendix [Cjl . 12 
As a consequence, the diagonal matrix element of the propagator also admits an TV-form. One virtue of the TV-form is 
that usually the expressions are much more compact. A drawback is that the functions I n ... rn do not directly apply 
for expressions written in TV-form. 

For the diagonal matrix elements of the propagator, through third order in the derivative expansion and in TV-form, 
one finds: 

(x\G(z)\x) h = J ' AP.h(p + iD)(N -2ip»N ll N - Ap^N^N - Ap^N^N 2 + TV^TV 
-2i Pli N^N v N - 2ip ll N vll N„N - 2ip tl N ll N uv N - 2ip ll N vv N ll N 
-g^AT^TV 2 - -ip^N^N 2 - -ip^N^N 2 - 2ip„NF llu N u N - -ip„NF^ u N 2 
+8ip^PvP a N^N„ a N 2 + SippPvpaN^NNaN + !6ip^ Pa N^N a N 2 + 8ip^ Pa N^ a N 3 
+8ip^ Pa N^N u N a N + 0(D 4 )). (4.14) 



The corresponding expression for the fourth order terms is given in Appendix [CJ In these expressions there are no 
ambiguities related to integration by parts in p^ or z and so the formulas are essentially unique. The only remaining 
freedom is to reorder the covariant derivatives. 

2. Trace of the propagator 

In order to obtain the trace of a generic function of K, one can use 

Tr h f(K) =[ ^f(z)Tr h G(z). (4.15) 

The expression of Tr/j G(z) can be obtained by starting from (x\G(z)\x)h (eq. (14.141) ) and using integration by parts 
and the trace cyclic property to obtain a simpler form. Due to the presence of the factor h(p + iD), the integration by 
parts (with respect to the covariant derivative) and the cyclic property do not act in the usual way for the expression 
in parenthesis. Instead one can use the identity 

^A(p)h(p + iD)B(p) = J ^h(p + l D)e^ Pf) -A(p)B(p). (4.16) 

Here A(p) and B(p) are arbitrary operators which may depend on p^ (but not on d?). ^ is [£> M , ] acting only on 
A(p). On the other hand, acts on the p^ dependence in A(p) and B(p). This identity is proven in Appendix iDl 

Of course, if one is working modulo 0(D n+1 ), and AB = 0(D n ), the operator e ldPDA can be dropped, and the cyclic 
property works as usual. 

However, the expression for TThG(z) is more easily obtained from the relation 

Tr,,G(z) = -^-Tr* log(z - K), (4.17) 

az 



12 We do not have a proof of this to all orders (but have little doubt that it is so). It has been verified through fourth order in the derivative 
expansion. 
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using the compact expression for Tr^ log(z — K) to be given below (eq. (|4.23[1 ). An explicit calculation to third order 
gives: 



, . d d xd d p 

' l = ' (27T) d 



% + iD) (a - ip^NN^Nu - MppN^uNv 



+ip^F^ w NNN v + ippF^NNvN + -ip^F^N^NN 

8 26 

+jjiptf>vPaNNNpN va + WipuPvPcNNN^Na + —ip^PaNN^NN^ 

-2ip^NN Vfl N u - 2ip lx NN v „N lx + 44ip^p u p a NN lt N v N a + 0(£ 4 ; 

Also the matrix elements of the propagator can be recovered from the logarithm by using 

6 



(x\G(z)\x) h 



SX(x) 



Tr h log(z-K), 



but in this case the relation eq. (|4.16[) is needed to extract the factor SX(x) in STih log (2; — K). 



(4.18) 



(4.19) 



The effective action in Chan's form 



As follows from eqs. (|4.15[> and (|4.18|) . for a generic function of K, Tihf(K) requires an integral over p^ and another 
over z. Nevertheless, the parametric integration over z can be obviated in the special case of the logarithm. Tr^ log K 
is just the effective action. 



1. Tr h log(z - K) 

As it turns out (verified through four derivatives) the diagonal matrix elements of the propagator can be written as 

d d p 



(x\G(z)\x) h 



(2tt)° 



h(p + iD)^^- + h{p)C(z] 
dz 



(4.20) 



Here M.(z) is a multiplicative operator that admits an A- form, and C(z) is traceless (a sum of commutators). There- 
fore, 



Tr h log(z - K) 



d( . , f d d xd d p 

— log (z-0 / u tr 
2m J (2ir) a 



h{p + iD) 



dM{Q 



d( 



(4.21) 



The term with C(z) has dropped from the expression. Next we integrate by parts in £, this transforms log(z — Q into 
\/{z — £)■ The integrand is assumed to be well behaved at infinity (in particular, the branch cut of the logarithm is 
no longer present). Hence, we can switch from the contour T, that includes the spectrum of K and excludes the pole 
at £ = z, to a contour excluding the spectrum of K and including the pole at £ = z. This produces 



Tr/j \og{z — K) 



d d x d d p 
{2n) d 



tr 



h{p + iD)M(z) 



(4.22) 



Now Tr/j log(z — K) is written in Chan's form, namely, in A-form and without parametric integration on £. Note 
that the dependence on z is inessential as z can be absorbed in X. 
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Explicitly, 



M(z) = -logJV +p„p v N fl N v 

112 2 

-iPuNNyFuv - -ip^N^NF^ - -ip^PvPaN^NNa + -ip^PvPaN^NaN 



1 7 28 

N aa N„ v N + 

n Pp,Pv 



-—p^PvN^N^Fua - -PpPuNN^NFva + yP^JV^JV 2 ^ 

-^-PvPuNaNN^Fva - ^-PnPvNN^N a F ua - -p^NN^F^ 
y y y 

~P*&> N t> NN <* F *<* - ^PnP„N^N a NF ua + -p^N^NF^ 

2 o 
—gPuPvN F n«F va 

8 2 

+ -^P^PvPccPpN ll N va Ni3N - ^Pfj.PvPaPpN^NNapN - kPuPvPaPpNuvNapN 

+^P»p v p a ppN tl N iy N a Np + 0(D 5 ). (4.23) 

(The isolated term — log A*" = log(iV _1 ) is still considered to be in iV-form.) 

The form of M(z) is not unique, due to the cyclic property and integration by parts with respect to the covariant 
derivative. 

That Chan's form exists is not trivial, in the sense that it holds for the logarithm but not for generic functions of K . 
Chan's form was introduced in [40j |. Extended to six derivatives in j47j . to curved space-time in (36j . and to fermions 
in |48| . It is quite remarkable that it also exists in /i-spaces (in particular, at finite temperature). This is more so as 
we are not allowed to use two important tools of the original derivation by Chan (40j , namely, momentum average 
and integration by parts with respect to p M . This is forbidden due to the presence of the function h(j>), which is 
arbitrary. It is noteworthy that, unlike the original Chan's formula, our expression does not depend on the space-time 
dimension. This property is also shared by the heat kernel. Another difference with Chan's result is that the are 
contracted only with covariant derivative indices and not with other p^ . 

2. Traced heat kernel 

To obtain the traced heat kernel, eq. (I3.28[) . from the effective action, eqs. (I4.22j) and (|4.23[) . one can use 



Tr h e K = I ^-e z Tr h 



dz „„ 1 



2iri z — K 



- /^e"Tr,log( = -A-) 



dz f d xd p 



2ni J (2ir) d 



h{p + iD)M(z) . (4.24) 



Now the shift z — » z — p 2 implies TV — » (z — X) 1 in A4(z), and e z — > e z e p . Hence, M(z) becomes p-independcnt 
and the integral over momenta reduces to obtaining the following /i-dependent operators 

(P* ' • • P„ n )h = (^) d/2 f^d h (P + Wy-^Pn ■ ■ ■ P*. • (4-25) 
The B n in section ImBl are obtained in this way. 



23 



3. Reduction to Chan's form 



In what follows, we explain how eq. (|4.23|) is obtained. First, let us see how Chan's derivation [40[ can be adapted 
to the present case. Using eq. (|4.1I) . 



(x\ log(z - K)\x) h 



d d p 

j^yi Hp) log (z - (d^ + i Pfl ) 

d d p 
{2n) d 

d d p 



X) 



h(p) [log (TV" 1 ) + log(l - (2ip^D fl + Dl)N) + C] 



C d d n 1 / \n 



(4.26) 



where C denote commutator terms, which will vanish upon use of the cyclic property of trace. To second order in the 
derivative expansion 



T r/l log(z-X) = 
Using the relations 



d d x d d p 
(2n) d 



tr 



/^(log^- 1 ) - 2ip^D^N - DlN + 2p^ Pu D^ND u N + Op 3 ))] . (4.27) 



the trace can be written as 
Tr h log(z - K) 



d^logiN- 1 ) = 2p^N, 

^dpPlogiN- 1 ) = S^N - 2 PlxPu N 2 , 

d d xd d p x 
{2iz) d 



(4.28) 



-tr 



h(p)(\o E {N- x ) - l D^log(N- 1 ) - -D^dldlXogiN- 1 ) 



-2 PllPv D^D v N 2 + 2p tl p„D ll ND y N + 0(D 3 ) 



d d x d d p 
(2n) d 

d d x d d p 
(27r) d 



tr 



tr 



h( P )e-* D9P (log(TV-i) +p, l p„N fl N v + 0(D 3 ) 
h{p + iD)(\og{N- 1 )+p^p y N tl N u + 0{D 3 ) 



(4.29) 



This expression has the desired Chan's form. 

In order to obtain the expression of Ai{z) to four derivatives it is not practical to apply the previous method since it 
is not sufficiently systematic. A possibility would be to simply write down all possible terms that could appear in Ai{z) 
to fourth order, with free coefficients, and expand everything in powers of D^, including h(p + iD) — » h{p)e~ lDd , 
using the cyclic property, to match the terms in eq. (|4.26[) . Assuming that the p^ can be only contracted with covariant 
derivatives (but not with other p^) the number of terms is finite (since iV _1 is not allowed). However, the number of 
possible terms is too large (and it is easy to miss some of them when trying to write down all of terms) . 

The method that we have followed is partially constructive and partially guessing. Let 



_ — * 1 / \ 71 

A(z) = \ogiN- 1 ) ((2*^ + Dl)N) + C, 



(4.30) 



where C are suitable commutator terms to be fixed. From previous formulas, 



Tr fe log(z-^) = J ±E±P t T[h(p)A(z)] 
d d xd d p 



(2ir)° 



h(p + iD)e lDd " A{z)e 



-iDd" 



Hence, we have to choose C, if possible, in such a way that the operator 

M(z) = e* D9P A(z)e^ D9F 



(4.31) 



(4.32) 
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is multiplicative and in iV-form. To see how this condition reflects on A(z), let us define two first-order variations, 
namely, 

6 D : D^-t iSa^, 

S p : p^ -> p p + 8a^ (4.33) 

where Sa^ is an arbitrary constant c- number (common to both variations). Clearly, the condition that M.(z) is 
multiplicative (and so with the covariant derivative operators in the form [D^, ]) is that 

5 D M{z) = 0. (4.34) 

Using eq. (|4.32[) . this requirement translates into the following condition on A(z): 

(S D - S p )A(z) = 0. (4.35) 

In turn, this is just the condition requiring that A(z) must depend only on the combination + ip^. This property 
is manifest in the symbol log(z — (D^ + ip^) 2 — X), but is not automatically preserved by the derivative expansion 
with formal use of the cyclic property (which is needed to have an iV-form) . So we have to choose the freedom implied 
by the cyclic property (i.e., the commutator terms C(z)) to fulfill eq. (|4.35|) . 

What we have done is to expand A(z) in eq. (|4.30|) , but allowing all possible cyclic permutations for each term, 
with free coefficients (this is the guess). Such coefficients are then partially fixed by the condition of reproducing 
log(2: — (D^ + ip^) 2 — X), modulo the cyclic property, and by the condition in eq. (|4.35j) . This condition is easily 
implemented by means of the rules 

(Sd - Sp)D ll = iSa^, (Sd ~ S p )p ll = -Sa^, (5d - S P )N = 2Sa t _ l p tJ ,N 2 . (4.36) 

The corresponding M(z) obtained from eq. (|4.32|) is multiplicative. It can be written in a manifestly multiplicative 
form by moving the to the right, forming covariant derivatives. The remaining freedom in the coefficients is used 
to obtain a simple form for A4(z). The guess chosen works at least to four derivatives, and very likely also to all 
orders. We conjecture that Chan's form for general h(p) can be extended to curved space-times as well. 



4- Quenched approximation 

The Polyakov loop in the formulas, or more generally the explicit iD^ in h(p + iD), is needed by consistency, but it 
is also a nuisance. Here we study the effect of setting this explicit iD^ to zero by hand in an expression. We call this 
the quenched version of the expression. The quenched results will be incorrect in general, but still one can consider 
whether this approximation can be done consistently. 

For any operator, the same derivation leading to eq. (|4.5I) can be repeated putting the explicit iD^ to the right: 

(x\f{D,X)\x) h = J ^L h (p + i D)f(D,X)= J ^Lf(D,X)h(p + iD). (4.37) 

After quenching, by setting iD^ to zero, the two expressions yield two different (incorrect) results. In fact, for a 
Hermitian operator like e , the unquenched matrix element respects hermiticity, but the two quenched expressions 
do not (rather they are hermitian conjugate of each other). On the other hand, inside the trace, the two quenched 
expressions do coincide (with each other, but not with the exact one containing h(p + iD)) 13 

Tr h , q (f(D,X)) = J ^0P t r[h(p)f{D,X)] = J ^ftr [f(D, X)h(p)] . (4.38) 

This relation provides a concrete choice of quenched version of the trace of an operator f(D,X). In general this will 
not coincide with first computing Tr^ (f(D, X)) and then quenching, since the latter does not commute with the cyclic 
property or integration by parts. 



Inside_the trace and the integral over p M , h(p)f(D, X) = h(p)e i9PD f(D + ip, X) = h(p - iD)f(D + ip, X) = f(D + ip, X)h(p - iD) = 
f{D,X)h(p). 
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Next, we study whether the quenched version of the traced heat kernel (as defined from eq. (|4.38[) using e K ) satisfies 
a consistency condition like that in eqs. p. 281) and (|3.29[) , namely, 

( x \e K \x) h , q = ^^^, Ti h , q (e K ) = J d d xtT(x\e K \x) h , q . (4.39) 

This kind of conditions, and similar ones for the effective action, can be derived from the corresponding relation for 
the propagator (which of course holds in the unquenched case too) 

<«<•> 

As a matter of fact, eq. (|4.40l) is correct, just because the expressions admit an iV-form and so depend on z and X 
in the form z — X. (Note that z — X appears only in N, and any given block JV^...^ can be transformed into N 
by integration by parts. So one needs to consider only the X and z variations on each block N at a time. The two 
variations give the same result thanks to the trace.) 

The quenched version of the effective action, as defined in eq. (|4.38[) using log K, does not admit a Chan's form. 
This can be seen from eq. (I4.14[) . (Note that eq. (|4.18[) has already applied the cyclic property and the integrand 
shown there is not unique.) After applying the quenching prescription in eq. (|4.14l) . one finds that the term N^^N 
cannot be expressed as a derivative with respect to z, modulo commutator terms. 

We stress once more that, in general, setting iD^ — in a traced quantity, written in different ways related by 
the cyclic property, yields different results. This follows from eq. (I4.16[) . At finite temperature this is also clear from 
eq. (|2.20[) : when using the cyclic property, commutation with Q produces time derivatives which are missed if Q is 
set to unity by hand. Also, relevant contributions can be missed by quenching. For instance, the first contribution to 
the induced charge density, obtained by taking a variation with respect to the potential Aq(x) in the effective action, 
comes from the Polyakov loop. Other contributions coming from Fy, v (x) contain more derivatives. 

All this implies that the quenched approximation is rather dangerous and may produce uncontrolled results. There- 
fore, quenching should cither be avoided altogether or at least, a careful evaluation of the contribution coming from 
the terms neglected should be done. 

V. SUMMARY AND CONCLUSIONS 

We have developed a new technique to deal with diagonal matrix elements of generic pseudo-differential operators. 
The method applies at finite temperature, or more generally, to /i-spaces, i.e., spaces with weighted integrals over 
the momentum of the loop. The approach is based on extending the method of covariant symbols to such spaces. 
This allows to carry out a manifestly gauge covariant and Lorentz covariant calculation throughout. We conjecture 
that the approach can be extended to curved space-time as well. In this case, the Polyakov loop of the Levi-Civita 
connection is expected to appear in parallel with the gauge connection. 

The new technique is appropriate to carry out covariant derivative expansions, so we have applied it to the heat 
kernel and to the effective action in Chan's form. For the heat kernel we present results for the diagonal matrix 
elements to three derivatives (the fourth order terms have also been obtained but are too bulky to be included). For 
the trace of the heat kernel we present results to four derivatives. We also present, to four derivatives, the expression 
of the effective action of a generic bosonic Klein-Gordon operator in Chan's form (i.e., prior to momentum integration) 
valid in /i-spaces. 

We have briefly touched the connection with the real time formulation of field theory at finite temperature. That 
formulation is appropriate to treat time dependent aspects, or even problems related to non-equilibrium physics. Such 
connection, not made in this work, is clearly worth pursuing. 

In this regard, we emphasize that the use of generic weights h(p) in this work is not intended as a device to describe 
time dependent situations or a connection to the real time formalism. 14 We merely observed that, although Hm{pq) 
is an even function of pa, all the Chan- like formulas can be written equally well for any h(p) without assuming any 
parity property or, more generally, any special dependence on p M . We find this formal property remarkable. However, 
as shown in eq. (14. 9p , the condition of gauge invariance does introduce a requirement of periodicity, which in practice 



Nevertheless, the possibility of transforming eq. fli.44t into eq. 1 12, 481 1 suggests that the use of weights h(p) defined on complex paths 
could have a bearing on this subject. Successful use of complex weights along complex paths in hadronic physics can be seen, for 
instance, in the spectral quark model introduced in [49|. 
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we only know how to fulfill for Euclidean space-times with the topology of a (possibly degenerated) torus. It could 
be that the formal property is just a mathematical nicety, or it could signal a deeper property of the formalism. This 
is not known at present. 

As emphasized in Section [TT] the derivative expansion implies a resummation of the dimensional (or large mass) 
expansion and in this sense we go considerably further that the standard approach for the heat kernel. The derivative 
expansion has the virtue of being gauge covariant at each order separately. In addition, higher order terms are 
increasingly ultraviolet convergent. So, for instance, anomaly saturating effective actions can be computed in closed 
form using this technique [34l . [38| . The derivative expansion in field theory is one of the few systematic tools to 
go beyond the perturbative regime H3, [HI, an d provides guidance in modeling of effective Lagrangians in exact 
renormalization group approaches |52j. The applicability of this technique extends to external field configurations 
which are slowly varying on space-time. In the case of external gauge fields, these have to be weak since they enter in 
the covariant derivative to preserve gauge invariance. The range of validity of the derivative expansion as compared 
to exact calculations in concrete profiles have been tested, most recently, in [HH . As expected the truncated expansion 
works better for quantum fields with short wavelengths, although even outside this regime it does a good job in some 
of the cases studied (53|. In general the derivative expansion is expected to be asymptotic, like the semiclassical 
expansion in quantum mechanics, with which it is closely related. Therefore a naive summation of higher orders 
would not provide a convergent result. To carry out the expansion to fourth order is rather standard, as the number 
of terms quickly increases beyond that order (35l. \at\. In addition, in four space-time dimensions, the calculation to 
four derivatives accounts for the ultraviolet divergent contributions. These are the ones leading to power counting 
renormalizable Lagrangian terms. 

Matrix elements of operators acting in the one-particle Hilbert space, (y|^4|x) , correspond to propagation along single 
lines from x to y in the Feynman graphs (in position space). General Feynman diagrams can be constructed joining 
these lines with the vertices of the theory under consideration. In the present work we have only studied diagonal 
matrix elements (x|A|a;). This is appropriate to produce local Lagrangians, under the derivative expansion. However, 
such restriction is certainly a limitation. For instance, as noted in Section III Al we can compute in this way the 
expectation value of the induced current (J M (a;)). By the same token, sum rules of the type (T(D n J(x))(D m J(x))) 
can be also obtained from the local calculation. On the other hand, the correlation between two currents in two 
different points (TJ fl (x)J^(y)) is not directly accessible. Therefore, an extension of the techniques discussed here to 
non diagonal matrix elements would be of interest. 

Upon completion of this work, we have learned that an equation equivalent to our eq. (|2.40p has been found 
independently in [HiJ in their study of CP violation in the standard model at finite temperature. 
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Appendix A: Commutator expansion 



If an expansion can be defined by means of a bookkeeping parameter the corresponding coefficients are well defined: 
fW = c n^ n an( i c n does not depend on how the expression is manipulated. Unfortunately, this is not the case 
for expansions based on counting the number of commutators. For instance, consider operators / = /(A, B) given by 
linear combinations of products of the basic operators A and B. No particular algebraic property is assumed for A 
and B (other than the associative property) . Let us grade the terms of the commutator expansion of / by the number 
of [A, ] they carry. This is ambiguous. For instance 

B 2 A = AB 2 - [A, B]B — B[A, B]. (Al) 

The expression as a whole is of zeroth order (this is the leading order). However, the concrete zeroth, first and second 
order components are different in the left and the right hand sides of the equation. 

To remedy this situation, the ambiguity can be removed by choosing a canonical form. A concrete choice comes 
from imposing the following prescriptions: i) In the canonical form the expression is written as a linear combination of 
products of blocks of the type, A or [A, ] n B with n = 0, 1,2, . . . (that is, B, [A, B], [A, [A, B}}, ...). ii) The blocks A 
are placed at the left. Further, the A's at the left count as order zero, and each block [A, ] n B counts as order n. The 
right hand side of eq. (|Alj) is written in canonical form: the zeroth order is AB 2 , the first order is —[A, B]B — B[A, B], 
and higher orders vanish. 

Let us now show that the canonical form just defined, as well as the corresponding grading of terms, can be derived 
from a bookkeeping parameter using labeled operators. Namely, by counting powers of A in 

/ = f(A, B)^f x = /(Ai + X(A — Ai),B). (A2) 

Here A\ represents A placed at the left (position 1 with respect to the blocks [A, ] n B). For instance 

B 2 A -» B 2 (A 1 + X(A - A x )) = AB 2 + \(B 2 A - AB 2 ) = AB 2 - X([A, B]B + B[A, B]). (A3) 

To proof eq. (|A2I) in general, first note that Ai — Ai + i is just [A, } placed at position i. E.g., (A2 — A^)B 2 = 
A 2 B 2 — A 3 B 2 = BAB — B 2 A = B[A, B]. Then, if a block A is located at position n, one can write 

A-Ax=A n -Ax = (A n - A„_i) + (A n _i - A n _ 2 ) + • • • + (A 2 - Ai). (A4) 

Therefore, A — A\ is a sum of commutators and A in eq. (IA2[) just counts the number of commutators [A, ]. 

This counting is unambiguous and extends trivially to the case of more operators, /(A, B,C, . . .) if terms are still 
graded by the number of [A, ]. It is worth noticing that things are more complicated for traced expressions, due to 
the cyclic property of the trace. (For instance, position "1" becomes ambiguous.) 

At the end of Section IIII A 21 it was noted that choosing to put the Polyakov loop at the left or at the right and 
then setting it to unity gives different results and breaks hermiticity of the heat kernel at finite temperature. This can 
be seen in the heat kernel coefficient Ai in eq. p. 261) since £2 is not zero for = 1. More generally, consider eq. (|2.43[) 
with Q — > 0. The covariant symbol K does not break the symmetry but the momentum derivatives can only be taken 
to the right due to the presence of /im(po) and this breaks the symmetry. Repeating the calculation with /im(po) 
placed at the right gives a different result, namely, the transposed of the previous one. The same conclusion can be 
obtained from ordinary symbols. This can be illustrated with a simple example. Consider the operator (Do + A) -1 . 
We consider the two expansions, with Do moved to the left or to the right: 

' ' = N L - N L (D - D , L )N L + ■ ■ ■ , (A5) 



D + X D ,l + X + D - D ,l 

here Nl = (Z?o,l + X)^ 1 and Do,l is Do at the left. If we set now Do.l ~> in Nl the result is 

N — N(D Q — Dox)N + 0(N 2 ) = N + [D , N]N + 0(N 2 ), N = X^ 1 . (A6) 

A similar calculation with Do r —> in Nr gives 

N - N{D Q - D 0tR )N + 0{N 2 ) = N - N[D , N] + 0{N 2 ). (A7) 

So the two prescriptions differ by [Dq,N 2 ] + 0(N 2 ). On the other hand, inside the trace, the two prescriptions do 
yield the same result (as stated in eq. (|4.38l0 . 
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Appendix B: Momentum integrals at finite temperature 

Let 

(?WW--PjO = {^r) d/2 J -^^h M {jp -Q)e- Tp ~p^ 1 p ll2 ---p tln . (Bl) 

These integrals are needed to obtain the heat kernel expansion coefficients at finite temperature. They are not 
normalized to unity. In particular 

(!>=$>. (B2) 
The basic result comes from distinguishing spatial from temporal degrees of freedom: 

11 

- ^r<5 ili2 ... i2n( .^_ )m ^ ro . (B3) 

Here the symbol ^ 1 i 2 ...i 2 „ represents the symmetric sum of the (2n— 1)!! products of n Kronecker deltas (each term 
with weight one). E.g. 

Sijki = SijSki + ^ikSji + SuSjk- (B4) 
Besides, we have introduced the auxiliary functions 

Vm = (47rr) 1 / 2 i ro r ro / 2 J ^ h M ( Po - Q) e^p™, m = 0, 1, 2, . . . (B5) 

These are related to the functions £ n of eq. (|3.3p through the relations 

m n 

<P™ = E *" +m 2 (n - m)/2 4 m £n = E (-0" +m 2- (n - m)/2 C„m ^m, (B6) 

n— m— 

where 

m n 

= E C 'nmHn{x), H n {x) = £ C„ m X™. (B7) 

n— m— 

As matrices c' = c _1T . 

In order to compute the heat kernel to four covariant derivatives, we need (p^p^ ' ' 'P/O f° r < n < 4. Using the 
previous formulas one obtains 



1 

W = 2^1/2^0^1, 



(p^PvPa) = ^_ 3/2 ^(^j/fco + <VAo + cW^oKi - 2^ ^ ^ q0 ^ 3 ) ' 

+^/3^0<^a0 + ^a^/iO^o) $2 + J^/iO^i/O^aO^oCl j 7 (B8) 



The formulas in this appendix plus the first eq. (|3.3p written as 

= (4^r) 1 /2(„i)«2-"/ 2 / ^° fc M (po - Q) e^H n (V^ P o), (B9) 

hold if the function /im(po) is replaced everywhere by a more general weight function, h(j>o). No special property of 
^Af(po) has been used. 



29 



Appendix C: Formulas 

Covariant symbol of K through fourth order in the derivative expansion: 

R = x - VilVll + ix^ae + v^ v dl - \x^dl 

-\v»F vailfi dldl% -F^dldl + -F^F^dldl + Lx^dldldldl + 0(D 5 ). (CI) 
Diagonal matrix elements of the propagator through third order in the derivative expansion, in X-form: 

f d d p ( 
{x\G{z)\x) h = J —^hip + iD)^ -2iI h2 p^X fl -Ah 

- (8/1,1,3 + *Ii,2,2)PpPvXnX v + 2I 1 . 1 . 2 X fi X fl - 2iI hia p„F^X v 
8 8 2 

- -^hflPytPvPoiFuva + -ih,ZPnPvPaF V na ~ -^h.^P^Fy^ + 8Ul,4PpPvP a X 

-% (-24/1,1,4 - 8/l : 2,3) P^PvPaX^Xya - i (-24/i, 1,4 - 16/l,2,3 ~ 8/1,3,2) PpPvPaX fj, v X a 

4 

-% (-48/i,i,i,4 - 32/i,i,2,3 - 16/l,l,3,2 - 16/l,2,l,3 - ^h^l^PiiPvPaX^XyXa - -ihfiPpXpw 

4 4 

--^ih,3PnX Vflv - -ih,zPnX vvll + i (-4/i,i, 3 - 2/1,2,2) p^X^X^ - kiI\j.$p il X v X ilv 

-4ih,i,3P l uX v X vtl + i (-4/1,1,3 - 2/1,2,2) PnX^ v X v + i (-4/1,1,3 - 2/1,2,2) p^X vtl X v 

+i (-4/i,i,3 - 2/1,2,2) p^X uv X^ + i (-8/1,1,1,3 - 4/i,i,2,2 - 4/1,2,1,2) p^X^XyXv 

+1 (-8/1,1,1,3 - 4/1,1,2,2) PnXvXpXv + 1 (-8/1,1,1,3 - 4/1,1,2,2) PpXvXvXp + Op 4 )) . (C2) 

Covariant symbol of G(z) in TV-form, through third order: 
G{z) = N + iN ll d*-2ip li N li N 

+N^N - ^N^dpP + 2p^N fl N v dP+p^NdP+Pf t N Vfl NdP 
+p lt NF ltu N8S - Ap^N^N - Ap^N^N 2 

+iN^N v 8l + l -iN^ v N8l + UN^Ndl + \iN^ v Ndl + ^iNF^Nd? 

-HNpF^NdZ - \iN^ a dldld a - 2ip^N^N vv N - 2ip^N llv N v N 

-2i Pfl N vll N v N - 2ip„N vv N l _ l N - ^ip^N^N 2 - ^ip^N^N 2 - ^ip^N^N 2 

-2ip^NF^N u N - y Plx NF vlxv N 2 + ip^N^d^ + ip^N^d^ 

+ip„N v ^N a dPd a + ^N^Ndldl + ^N^Ndldl + \i P »N ua ^Ndld a 

+iPnNF^N a dPd p a + ^ip^NF^ a Ndld p a + ip^F^NdPdP - Up^p v N^N v N a d a 
-2i PllPy N^N va Nd a - 2ip llPv N li N av Nd* - Up^p v N^NN a d a - Aip^p v N^N a Nd a 
-2ip^p u N m N u Nd a - 2i Pll p v N all N v Nd a - ^ip^N^ a N 2 dr - ^ Pu N^N 2 d a 

~ip^N a ^N 2 d a - 2i Pll p v NF^N v Nd p a - ^ip^NF^ a N 2 d a - 2tp^p v N ll F va N 2 d p a 
-2i Pl ,p v N^NF ua Nd p a + Sip^PvP a N^N u N a N + Sip^p vPa N ^N va N 2 

+8ip»p„p a N tiV NN a N + mip^ Pa N^N a N 2 + Sip^PaN^N 3 + 0{D 4 ). (C3) 



30 



Fourth order of the diagonal matrix element of G(z) in TV-form: 
(x\G(z)\x) hA = f ^P h (p + iD)[ 



d d p 

■nyp - 

2 „ „ „ 2 „ „ „ 2 „ „ „ 1 „ ~ 1 



N^NuuN + -N^ V N V N + -N^N V N + -N^N^N + -N^N 2 + -N^N* 
+\n iivvix N 2 + ^NF MMV N V N + 2N M F fu ,N„N + ^N^F^N 2 + ^NF^F^N 2 

8 8 8 

-^PuN^NuaaN 2 - -p^PvN^N^N 2 - -PvPuN^NvavN 2 - A PliPv N^NN aa N 

-8p^N^N a N a N ~ 8p^N^N aa N 2 - 4 PflPv N^ a N v N a N - 4p tl p v N lux N a N 1/ N 
-Ap^N^N^N 2 - ipuPuN^N^N 2 - 4 PfJlPv N ati N v N a N - i PtlPu N al ,N a N v N 
-A Pli p„N ai ,N va N 2 - A Pll p v N a ^N alJ N 2 - ip^N^N^N - Ap^N^N^N 2 

8 16 o 8 16 9 

-TpuPuN^aNNaN - —p^N ^ ua N a N 2 - -p^p v N^ av NN a N - —p^N^N^ 2 

8 16 8 16 

-^PuN^ aa NN v N - —p^PuN^NvN 2 ~ -p„p„N a ^NN a N - —p^p v N allv N a N 2 

8 16 8 16 

—£ptf>uN aiM NN v N - —PuPuN^NvN 2 - -p„p„N aafl NN u N - —p„p v N aaii N u N 2 

^p^N^aN 3 - 2p^p iy N fiaua N 3 - 2p il p v N iiaav N 3 - 2 Pll p v N a ^ a N 3 
-2p^p u N aliau N 3 - 2p ti p u N a N 3 - Ap^ L p u NF m N v N a N - 4p^p v NF^ a N a N v N 

8 16 

-4p„p v NF m N„ a N 2 - 4p^ Pu NF^ a N au N 2 - -p^p v NF^ a NN a N ~ —p^NF^N^ 2 

~P»p v NF aflot NN v N - ^P^NF^N.N 2 - p^p u NF^ iava N 3 - p^NF a ^ a N 3 
-AppPvNpFvaNNaN - %P^p v N^F ua N a N 2 - ^p^p v N^F ava N 3 ~ 4p^p y N li NF ua N a N 

~^P^PuN^NF ava N 2 - 2p fl p y N F^ a F ua N 3 + l§p^p v p a ppN^N u N a N p N 

+l6p li p„p a ppN ll N v N af} N 2 + IQp^pappN^N^NNpN + 32p llP „p aPfj N t ,N l , a NpN 2 
+16p»PvP*Pp%N uaf) N 3 + lQp^p v p aP0 N^NN a N p N + lGp^p„p a ppN^NN aP N 2 
+Z2p li p v p a ppN liV N a NN f} N + 64p^p aPp N^N a N f3 N 2 + 48puP v p a p p N ^N af} N 3 
+lQpnPvP a Pi3N^ a N 2 NpN + 32 Pf , PlJ p a p N^ a NNpN 2 + 48p^p a ppN^ a N p N 3 
+l6p^p a PpN^ a pN 4 ) . (C4) 



Appendix D: The cyclic property in /t-spaces 



In order to prove eq. (I4.16[) . we can assume, without loss of generality, that A{p) — Aa(p) and B(p) = Bb(p), where 
the operators A and B do not depend on p^ and a(p) and b(p) are c-numbers (i.e., they commute with everything, 
except d?). 
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/ J^S A (p)*(P + i D ) B (p) = J ^Ah( P + l D)Ba(p)b(p) 

= J ^Ah(p)e-* 9PD Ba(p)b(p) 

= J ^Hp)Ae- ldPD Ba(p)b(p) 

= J ^h(p + l D)e WP ^ABa(p)b(p) 

= J ^I h (P + ^)e WP£>A A(p)B(p). (Dl) 
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